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t^ 1 ■ Abstract 

■ Present quantum Monte Carlo codes use statistical techniques adapted 

to find the amplitude of a quantum system or the associated eigenvalues. 
, Thus, they do not use a true physical random source, ft is demonstrated 

that, in fact, quantum probability admits a description based on a spe- 
<^^> cific class of random process at least for the single particle case. Then a 

first principle Monte Carlo code that exactly simulates quantum dynamics 
can be constructed. The subtle question concerning how to map random 
choices in amplitude interferences is explained. Possible advantages of 
this code in simulating partial histograms observed in particle diffraction 
' experiments are discussed. 
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Predictions in quantum mechanics of pure states come from amplitude cal- 
culations. Despite the fact that information concerning particle dynamics is 
stored in amplitudes, no experiment can be performed which directly measures 
the value of an amplitude. For this reason, probabilities have a rather unusual 
source in quantum theory for they come from squared amplitudes rather than 
from a random source. In this sense, they must be considered as "a posteriori" 
probabilities, since they are obtained after the physical problem has been solved 
by using amplitude calculations. In contrast, in classical physics particle dynam- 
ics provided by stochastic processes are calculated by using " a priori" probabili- 
ties, obtained from specific physical models. Thus, Monte Carlo methods, when 
applicable to quantum theory, are conceptually different from those used to 
simulate classical dynamics. In fact, two classes of Monte Carlo methods, suit- 
able to meet the requirements of the quantum theory, were devised. One type 
uses Monte Carlo integration techniques to solve numerically the Schroedinger 
equation pQ. This class provides the correct wave function for complex many- 
body problems using some available mathematical techniques [2], [3], [I]. The 
other class uses calculated quantum probabilities to simulate particle motion 
subject to quasi-classical scattering [5], [6]. In this case the quantum problem of 
a basic process is analytically solved and the probabilities obtained are used to 
simulate complex situations involving an ensemble of these basic processes. In 
both classes, the primary goal is to obtain some amplitude prior to getting or 
to using probabilities. Consequently existing quantum Monte Carlo codes do 
not simulate quantum dynamics using a priori probabilities as is usual for real 
stochastic processes. Both theories make use of probabilities but they generate 
completely different Monte Carlo codes. This reflects the intrinsic reality-free 
interpretation of quantum theory which denies any use of sets of classical tra- 
jectories (like those provided by classical Monte Carlo simulations) as a reliable 
description of quantum phenomena. 

In this sense no devised stochastic dynamics has been able to fully simulate 
quantum dynamics. Notwithstanding direct calculation of quantum probabili- 
ties expressed as a sum of histories forming a stochastic-like process was per- 
formed by Marinov [7] . He proved that squared amplitudes can be written as the 
path integral involving a product of transition matrices defining histories in a 
classical-like phase space. Marginal calculations, obtained from partial integra- 
tion over momentum variable, lead to histories which are written as a functional 
of classical trajectories described by a Langevin equation. His formalism mimics 
completely a (classical) stochastic process, except that resulting transitions ma- 
trices are quasi-probabilities instead a positive-defined transition matrix. The 
real-valued character of his transition matrices provides some negative- valued 
paths that are necessary to produce quantum interferences. This stochastic-like 
process cannot be directly modelled by a Monte Carlo code since the paths 
are constructed from quasiprobabilities which cannot be used to draft random 
numbers, the seed of a first-principle-defined Monte Carlo code. This is a clear 
evidence that a simple stochastic process cannot model any quantum dynam- 
ics, as already demonstrated by Baublitz |8J. In contrast to this immediate 
interpretation of Marinov 's results a recent work of Skorobogatov and Svertilov 
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[5] has given support to the idea of a priori probabilities generating quantum 
processes. They demonstrated that quantum probabilities for an isolated two 
level system correspond to a particular kind of Chapmann-Kolmogorov equa- 
tion presenting both a non-Markovian character and a specific discrete jumping 
process. The formal stochastic process they obtained was entirely deduced from 
quantum dynamics with no additional hypothesis included. Thus, it appears 
that is possible to mimic some types of quantum problems by a suitable class of 
formal stochastic process. These results rise some relevant questions that ask 
for a solution. It is necessary to solve the apparently contradictory results of 
Baublitz and Skorobogatov and Svertilov, in order to understand which class of 
quantum problems admits a formal stochastic representation and to determine 
whether is possible or not to write a quantum Monte Carlo code based on a 
priori probabilities. Positive answers to these questions would probably lead 
to new possibilities on computation techniques of quantum problems and, at a 
more fundamental level, could give some clues about the still unsolved questions 
concerning the subquantum world. 

In this work, we touch upon and partially answer these questions. We present 
a specific class of formal stochastic process whose time evolution forms a chain 
mediated by quasiprobability matrices. Then it is demonstrated that Mari- 
nov functional is a particular instance of these matrices. In sequence, a priori 
transition probabilities that generate these real-valued transition matrices are 
constructed and an associated Monte Carlo code can then be written using 
them. Differently from usual quantum Monte Carlo calculations the obtained 
code is fully based on a classical algorithm without the need of make refer- 
ence to amplitude calculations although still capable of fully simulate quantum 
dynamics. In the code described here, quantum probabilities are obtained as 
an excess probability above a proper level, to be defined below. The delicate 
question of how to explain amplitude interference in terms of a proper choice of 
random numbers is discussed and its possibility demonstrated for the first time. 
Since a Monte Carlo code is the ultimate and definitive description of a random 
process, the existence of such a code, capable of model quantum processes, puts 
in order the need of a stochastic justification to quantum theory 

This work meets the requirements on presenting a first principle Monte Carlo 
code based on a formal stochastic dynamics. It has, in this sense, the operational 
goal of solve problems in quantum mechanics, but also pursues the challenge of 
to give clues to important questions involving the very nature of the subquantum 
world. In what follows we develop a model of specific stochastic processes pre- 
senting dynamics mediated by real-valued transition matrices. Then we present 
arguments that proves how a particular form of Marinov's results match this 
model. As a result a positive-defined transition matrix is derived from quantum 
theory and a complete chain Monte Carlo process is then trivially obtained. 

A Monte Carlo code works upon choices of random numbers, forming histo- 
ries. A process is simulated by the use of specific rules defining the way choices 
are selected as valid and on the way histograms are calculated. Thus the bare 
result of a Monte Carlo code is a set of valid counts. The difficulty in treating 
quantum events using a Monte Carlo code comes from the possibility of interfer- 
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cncc effects. For consider two quantum states labeled .4 and B, with amplitudes 



* (A) and * (B) and let p = |* (A) + * (B)| 2 the total probability for the state 
n = \& (A) + This probability can be null even if the probability for 

individual states, (A)\ 2 and (B)\ 2 , is not null. This means that we cannot 
take the number of counts generated by event described by state A alone and 
separately sum to the number of counts generated by event B in order to get 
the histogram for the total probability. 

In order to write a code adapted to quantum mechanics it is necessary to 
know how to handle, generically, process of type = s + (— s), where s is 
a number of counts in a histogram. Since there is no "negative" histories it 
is immediate that no code, based on direct choice of random numbers, may 
generate this result. In what follows we introduce the main argument of this 
work. It proves that, in an extended probability space (EPS), a code providing 
this result in fact exist. Let's consider the case s = kw, where both k,w G 
[— 1, 1]. Choose numbers g G [0, 2] and X,v G [—1, 1] such that k = 1 — g and 
w = A — v. Define K = (1 — g) I, where I is the two-dimensional identity matrix 
and write a set of vectors, W, Po and V, related by 
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The product giving s is equal to (KW) . Since the elements in Po and V are 
positive and sum to one they can be treated as probabilities. In similar way we 
define a quasiprobability vector S from the number s and get, from vector W, 
the result: 

S= ( 3 ) =KW = KPo~V + g ( , V , V )Po 
\—s J \1 — vl— vj 

^P^S + V=( 1 -^ 1 - V) ° V )po 

The columns in the last matrix sum to one. Their elements are positive if 
g(l — v) < 1 and gv < 1. Both relations demand that g < 2, a condi- 
tion naturally fulfilled. For K > (g < 1) any value of v G [0,1] is al- 
lowed. For K < (1 < g < 2) allowed values of v are restricted to the interval 
((1 — g) g^ 1 , g^ 1 ] ■ Within these restrictions, a positive matrix M is obtained 
that satisfies P = MPo. Therefore equation P — MPo maps directly into a 
formal stochastic process describing a system possessing two states subject to 
random transitions between them [lOj . In this process the system in the state 
indexed as 1 and possessing probability Po\ may transit, with probability gv, 
to the state indexed as 0; if state is the system may transit, with probability 
g (1 — v), to state 1. Total probability for state after this random process is: 
Po = (1 — g (1 — v)) Poq + gvPo\ = (MPo) . A first principle Monte Carlo 
code, using these transition probabilities and standard programming technics, 
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may model vector P. Now run H histories of this code. The number of them 
ending on state is P$H = (s + v)H and on state 1 is P\H = (1 — v — s) H . 
Assume s > 0. In discarding the first vH histories ending in state results in 
the value Hs. This is a permitted operation in a Monte Carlo code for it is 
one type of criterion used for validation of a history. Thus the number s > 
can be properly simulated using a first principle Monte Carlo code. If s < its 
absolute value is found from similar reasoning applied to state 1, from which 
(1 — v) H histories are discarded. 

We now treat the case = s + (— s). Modelling the event s + (— s) = 
kw + k(— w) means that a random choice must select event s or event — s and 
run a history for the chosen process. The resulting histogram is just addition of 
the obtained histograms. Select the same value v for both processes and consider 
a set of H histories. Select event s or — s with equal probability. The histogram 
modelling event s equals to (S + V)H/2 and event —s has a histogram equal 
to (—S + V)H/2. Total histogram, which describes the occurrence of both 
events, equals to (S - S + 2V) H/2 = VH. Thus, in discarding the first (VH) 
histories results in exactly null counts of the histogram to state 0. That is, this 
code can model the algebraic sum = s + (— s), in the sense that well-defined 
algorithms are associated individually to these events and its joint occurrence. 
Using a proper criterion for validation of histories, a result of null counts is 
generated, we now show how to use this method in a form appropriate to model 
real quantum problems. 

The starting point is the time-slice Feynman formula 

— r / ■ \ 

* (a; ' t] = (^fe) T J exp u An ) *° {xo) n dxi 

J \ / l = Q 

which is know to converge to the Schroedinger wave function when t/n = e — ► 
. In this equation the symbol A n is the discrete-time action written as 
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with end point x n — x. The point chosen to calculate the potential rises subtle 
questions [TT] and the simplest version, but not necessarily the more precise 
one, is assumed here. The squared amplitude, written in non-dimensional units 

x — » xy^he/m: U — > h/eU, becomes 

I* {x)S - f J exp (i [A n { Xl } - A n {*{}]) * (ar ) % (x' Q ) [] dxM 

-00 1=0 

As before it is convenient to recast this equation by setting xi — ui + .5wi and 
x[ = ui — .5wi. Expansion of the argument of the exponential on these new 
variables puts the action in a new form: instead the two-time chain present 
in the kinetic energy term, a three-time chain on the variable u appears in 
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association with a single-time dependence on w. The result is a odd series in 
this variable of the form 



£ 

1=1 



(ui+i + ui-i - 2u{) wi + U [ui + ^wi J -U (ui- ^wi 



I Hi </,,! - U ( U + T^" '- j + / ' (u Q - -W 



Partial integration, for indices greater than zero and at independent time slices, 
of the exponential that appears in eqn([T|) is admissible. Due to the odd parity 
in w, it results in a real-valued kernel for the probability, valid for individual 
time slices, and given by 



Ki [ ui + i + - 2ui + —U (ui) ; u t 



cos 



ui + i + ui-i - 2ui + —U (ui) ) wi + 2 



k 



W 



2k+l d 2k+l U{ui y 



J(2fe + l)l 



du* k+1 



dwi 



In the above equation we have a kernel valid for short time intervals. Time 
slice zero deserves a special issue. The first coefficient of wo is just the velocity 
written in terms of variables u. Since there is no chain term in wo, integration 
in this variable can be performed separately and has the form 

oo 

-!- j e i(«i-«o)«$ ^ Uf) + e-i^uo+H** (u + ^ e lV ( Ua ~? w )dw = W (« 0l ui - u ) 



This leads to a modified Wigner function at time zero: an additional phase 
term changes the value of the initial wave function. Its existence does not make 
invalid the line of reasoning presented here so we postpone discussion on its 
meaning. Putting together the last results it is clear that quantum probability 
can be written as a path integral of real valued kernels possessing as initial 
condition a Wigner function . Explicitly we have 



°° 71 — 1 

\^{x)J 2 dx = dx j JJjfftti+i+uj-i 



2u, 



d 

dm 



\ " 1 

-U (ut) ; ui j W (u , wi) Y\ du i 
1=0 

which is a particular case of a more general result obtained by Marinov [7] who 

get similar expression but using integrals in the phase space instead. We now 

approximate these integrals by a discrete sum on a finite spatial grid. In this case 

r n-i I 
Aa; II KiAui 
1=2 



each path consists of a product of type 



WAuoAm = k ef W 
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where fe e / = Ax Y[ KiAui and W = WAuqAui. The sum over all paths 

1=2 

becomes a sum over all values of the product k e fW . Both are real numbers 
with absolute values smaller than one. Thus the formalism of extend probability 
space just described is applicable to every path entering in the above sum. A 
Monte Carlo code that selects a path by chance and simulates the product 
k e fW in the extended space generates a count that approximates to the result 
of eqn© . 

In fact more information is available by eqn(??) defining the kernels. Since 
the kernel k e f is factorable, the first product in a path has the form K\W = 
[Mi (Po — V)] = (Pi — V) where the last equality comes from the fact that 
the reference vector V (of a given path) is invariant under application of the 
matrix M. The second product has the form K2K1W = K 2 {M\Pq — V) = 
(M>}M\Pq — V) . Recursive application of this method on the product involving 
the whole chain of a path leads to 



\*(x) n \ 2 dxf*(P n -V path ) = 



Jl M t (P - V path ) 



1=2 
.path 



(3) 



This result clearly express that probabilities in quantum theory are, in the EPS 
formalism, a difference of probabilities. Using the invariance property of V , we 
get P„ = Yi MiPq. This equation represents a classical probability process in 

path 

the EPS. Quantum effects comes from two effects: a renormalization of the prob- 
ability and from the swapping of states in the EPS. Notice the argument in the 

transition matrices M (ui+i + — 2ui + j^U (ui) , uij is a time-difference 
equation 

d 

ui+i + - 2ui + —U (ui) = yi 
dm 

It is possible to interpret this expression as a Langevin equation possessing 
noise source y (with uniform probability distribution) and linking the various 
time slices, in such a way that the transition matrix is function of a stochastic 
force through the random variable y. This way, total probability at time t is 
given by 

p (' T )= E U M (yi^i)P 

{path} I 

where all paths must finish at x. Thus, in the EPS formalism, a particle in 
quantum motion follows a simultaneous two-step process. It is subject to a 
classical random process displayed in the Langevin equation. The same random 
variable that drives the stochastic force performs the swapping process in the 
EPS. Initial probability is calculated from the addition of a proper reference 
value to the Wigner function that makes the result everywhere positive. The 
associate Monte Carlo method, with H histories, has the following structure: 
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• a spatiotemporal grid is mounted and approximate Marinov kernels cal- 



• from the knowledge of these kernels an appropriate reference value is de- 
termined for each path; 

• transition matrices are calculated; 

• initial joint probability P Q is constructed from the Wigner function; 

• initial positions at time zero and one are selected from the initial joint 
probability Pq; 

• a random number is selected; 

• it determines the next position, calculated form the Langevin equation; 

• in the same time step swapping of states are calculated using the proba- 
bility values displayed in the matrix M (y; u); 

• sequential choices of random numbers lead to a path, all them restricted 
to terminate at position x; 

• normalized histogram for the states, at each ending point, is generated 
from the counts; 

• quantum probability is get from the excess probability above the appro- 
priate reference level; one state always has probability above its reference 
level and the other one below; 

• another end-point is chosen. 

Thus, as anticipated, a Monte Carlo code capable of calculate quantum 
probabilities and using first principle transition matrices can be constructed. In 
addition, the general form of the short-time Feynmam propagator led to a naive 
interpretation for the associated quantum probability path integral, by which 
space dynamics (conveniently described by a Langevin equation) matches the 
swapping process involving the states of the EPS, both process driven by one 
single random variable. 

Due to the dependence of the eqn(??) on the value of the potential energy 
at time zero, the exact definition of the initial quasiprobability rises the still 
elaborate question on (quantum) state preparation [T2] , since this procedure also 
depends on the way a potential changes the phase of the initial wave function. 
This technical question affects the exact value of the zero-time amplitude but 
does not introduce conceptual changes in the algorithm proposed here. Another 
important feature is the non-local character of the expression for the quantum 
transition quasiprobability, the kernel K, expressed by the presence of high- 
order derivatives in its argument: 



culated; 
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and this imply that information concerning the existence of a potential is spread 
over the whole space even for very well localized potential profiles. Then it might 
be possible to capture some of the non-local effects existing in quantum theory 
as in the very relevant case of EPR-Bell experiments. Furthermore note that 
eqn((3]) is valid even for time-dependent potentials which just adds a index to the 
potential energy calculated in the spatiotemporal grid but does not introduces 
additional chaining effects. Thus the equivalence of the code proposed here and 
single particle quantum dynamics in one dimension is complete for a wide class of 
problems, including time-dependent potentials and arbitrary initial conditions, 
possibly with inclusion of technical complications induced by state preparation. 

The scenario is clear. Likewise the explanation given for a stochastic case, 
in the present interpretation of the results of Marinov, quantum mechanics re- 
sembles a ensemble of systems existing in a state of the EPS at values consistent 
with the reference level (the "vacuum"). Measurements makes sense only for dif- 
ference of probabilities, relative to vacuum values, as shown in eqn([3]). Here, the 
interpretation of quantum probabilities gets more involved than the usual one. 
Quantum probability is a positive-definite quasiprobability which describes the 
likelihood of a specific internal state (the observable one), at some place. This 
means, in the present interpretation, that quantum theory cannot give com- 
plete account of the actual physical state of a particle. In fact, the exchange 
of internal states during the time evolution of the stochastic dynamics leads to 
the effect of interferences of quasiprobabilities provided by the amplitude for- 
malism. Only one internal component of the stochastic process is observable 
because the other has probabilities necessarily below the associated vacuum 
level, thus always presenting a negative quasiprobability, so the particle appear 
to present no internal state at all. This effect hidden the stochastic nature of 
the process, when described by using amplitude calculations, because the action 
of individual transition probabilities no longer exists in the diagonal evolution 
of the quasiprobability vectors. In consequence the formalism of quantum the- 
ory have less information about the actual dynamics of the particle than the 
stochastic one. It may be argued that no measurements can detect this effect, 
but as will shown bellow, for a specific class of experiments, it seems that the 
stochastic model may capture more details than the amplitude formalism. This 
way we have here a clear and precise description of the importance of trajec- 
tories in quantum theory, a rather involved question on the very foundations 
of quantum mechanics. When simulated by a first principle Monte Carlo code, 
quantum trajectories appear to be real in the same sense they are, for instance, 
in codes modelling diffusion. All paths are tested and included in the statistics. 
All paths starts from a initially distributed region of the space and must get 
the same final observation point. In this aspect we have a complete equivalency 
to classical codes, describing by usual stochastic dynamics. The vacuum level, 
which determines measurements results, matters in the quantum case and has 
no meaning in classical dynamics. Besides, much more trajectories are used 
to get the probability at a point, namely those involving the extra (unobserv- 
able) degree of freedom. These two elements, which are correctly considered in 
the logical structure of the extended probability space theory, make quantum 
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trajectories to appear as possessing no physical reality. That's the essence of 
quantum interferences, in the present formalism. 

As far as the arguments presented here go, we have demonstrated that non- 
relativistic quantum mechanics of one particle admits the existence of a formal 
Monte Carlo code based on a priori probabilities. In common to other Monte 
Carlo applications, this code have full advantage over analytical methods when 
complex mechanical problems arc treated. No physical justification was given to 
the random choice based on extra degrees of freedom nor to the mathematical 
properties of transition probabilities either. Despite this limitation, it permits 
a very general modeling of quantum phenomena, and represents a new class 
of quantum Monte Carlo method. It may be especially useful in modelling 
single hits experiments like that performed by Tonomura, Endo, Matsuda and 
Kawasaki |13) . They shot single electrons past a double slit apparatus and ob- 
served the formation of the interfering pattern. In this case squared amplitudes 
give only the final histogram valid for a large number of hits. The code proposed 
here, when extended to the two dimensional case, can explain classes of inter- 
mediate histograms, valid for arbitrary number of detected particles, indicating 
that, in this case, more information is get using a stochastic process than using 
amplitude calculations. 

The nature of the noise source deserves a special issue. This question ex- 
trapolates the model developed, which is consistent with information contained 
uniquely inside the amplitude formalism. Its physical reality relies on the strong 
sensitivity of Marinov quasiprobability on the potential energy as well as the 
possible existence of extra degree of freedom, both of which demands proper 
justification on experimental grounds. Another point touches the way simu- 
lated trajectories are interpreted. It is well known that trajectories in quantum 
mechanics are continuous, as shown by Feynman [T7] and the same does hap- 
pen in ordinary Brownian motion, this fact granted by the Lindeberg condition 
[15] . Despite these similarities quantum dynamics cannot be explainable by any 
stochastic force associated to simple classical noise sources unless some nonclas- 
sical procedure are used such as a negative diffusion coefficient [TB] . The extra 
degrees of freedom introduced here and vacuum renormalization of the stochas- 
tic probability, which are acceptable logical elements inside a classical reasoning, 
appear to complement the properties of pure classical noise sources in order to 
reproduce quantum effects. This way we have got a consistent conciliation of 
both models and explained the origin of its widely different behavior as well the 
common probabilistic nature, namely the dynamics provided by the Langevin 
equation. This argument effectively solves the apparently contradictory results 
presented in ref. [8] and [9] . Since it is now possible to interpret quantum theory 
through a proper stochastic view the possibility of a real noise source driving 
quantum phenomena cannot be discarded at all. Such possibility demands addi- 
tional studies in more complex situations, as those involving the Dirac equation 
or many-particle systems, in order to explore more deeply the possible existence 
of a hidden stochastic mechanism in these cases as well. Probably, additional 
properties of the transition matrices may appear when these cases are considered 
giving new formal clues concerning the nature of the vacuum source underlying 
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quantum phenomena. 
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